Astronomy & Astrophysics manuscript no. paper 


© ESO 2008 


October 6, 2008 





Large-scale dynamos in turbulent convection with shear 

P. J. Kapyla 1 , M. J. Korpi 1 , and A. Brandenburg 2 

1 Observatory, Tahtitorninmaki (PO Box 14), FI-00014 University of Helsinki, Finland 

2 NORDITA, Roslagstullsbacken 23, SE-10691 Stockholm, Sweden 

Received 2 June 2008 / Accepted 30 September 2008 

ABSTRACT 

Aims. To study the existence of large-scale convective dynamos under the influence of shear and rotation. 
Methods. Three-dimensional numerical simulations of penetrative compressible convection with uniform horizontal shear 
are used to study dynamo action and the generation of large-scale magnetic fields. We consider cases where the magnetic 
Reynolds number is either marginal or moderately supercritical with respect to small-scale dynamo action in the absence 
of shear and rotation. Our magnetic Reynolds number is based on the wavenumber of the depth of the convectively 
unstable layer. The effects of magnetic helicity fluxes are studied by comparing results for the magnetic field with open 
and closed boundaries. 

Results. Without shear no large-scale dynamos are found even if the ingredients necessary for the a-effect (rotation and 
stratification) are present in the system. When uniform horizontal shear is added, a large-scale magnetic field develops, 
provided the boundaries are open. In this case the mean magnetic field contains a significant fraction of the total 
field. For those runs where the magnetic Reynolds number is between 60 and 250, an additional small-scale dynamo 
is expected to be excited, but the field distribution is found to be similar to cases with smaller magnetic Reynolds 
number where the small-scale dynamo is not excited. In the case of closed (perfectly conducting) boundaries, magnetic 
helicity fluxes are suppressed and no large-scale fields are found. Similarly, poor large-scale field development is seen 
when vertical shear is used in combination with periodic boundary conditions in the horizontal directions. If, however, 
open (normal-field) boundary conditions are used in the a;-direction, a large-scale field develops. These results support 
the notion that shear not only helps to generate the field, but it also plays a crucial role in driving magnetic helicity 
fluxes out of the system along the isocontours of shear, thereby allowing efficient dynamo action. 
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1. Introduction 

It is evident that the Sun possesses a large-scale mag- 
netic field. The most likely origin of this field is a hy- 
dromagnetic dynamo working within or somewhat below 
the turbulent convection zone (e.g. Ossendrijver 2003). The 
solar dynamo is generally thought to rely on large-scale 
shear to produce toroidal magnetic field from poloidal field 
and a process called the a-effect which produces poloidal 
field from toroidal field. In the turbulent mean-field dy- 
namo picture an a-effect can arise due to helical motions 
in a stratified fluid in the presence of either rotation or 
shear, or both (Rogachevskii & Kleeorin 120031 Rudiger & 
Kitchatinov 2006; Radler & Stepanov 12006!) . Other candi- 
dates for regenerating large-scale magnetic fields in shear- 
ing turbulence include the incoherent afi-effect (Vishniac 
& Brandenburg 119971 Proctor |2007[) and the shear-current 
effect (Rogachevskii & Kleeorin |2003| \MB$ ■ 

Numerical simulations of convection have been used to 
study dynamo action in local (e.g. Cattaneo & Hughes 
|2TjTJ61 Tobias et al -120051) and global (e.g. Brun et al. 120031 
Browning et al. 120061 Brown et al. I2007P settings. So far, 
however, with the exception of the last two papers, these 
models have had a hard time generating appreciable large- 
scale magnetic fields, although all the ingredients such as 
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rotation and stratification necessary for an a-effect to occur 
have been present. Possible reasons for the lack of large- 
scale dynamo action in the numerical models is arguably 
the lack of scale separation between the energy carrying 
scale and the scale of the domain and the absence of large- 
scale shear in the system. In stratified convection, the lack 
of sufficient scale separation is difficult to tackle due to the 
prohibitive computational challenges involved when trying 
to cover many more vertical scale heights in the domain. 
Adding shear, however, is easier to accomplish, which is 
one of the principal aims of the present paper. 

An additional point concerning large-scale dynamo ac- 
tion in turbulent convection is the possibility of catas- 
trophic quenching due to small-scale magnetic fields. It has 
long been argued that small-scale fields cause catastrophic 
quenching of the large-scale dynamo (e.g. Vainshtein & 
Cattaneo I1992p . This sort of quenching, however, only oc- 
curs in specific circumstances, i.e. when there is no mag- 
netic helicity flux out of the domain. This behaviour can be 
understood in terms of magnetic helicity (Brandenburg & 
Subramanian 2005a) which is a conserved quantity in ideal 
MHD. If magnetic diffusivity is finite, the magnetic helicity 
can only change diffusively under the aforementioned spe- 
cial circumstances that lead to resistively slow saturation 
of the large-scale field. 

If, however, magnetic helicity is driven out of the do- 
main, catastrophic quenching can be alleviated and the 
large-scale dynamo is expected to saturate near equipar- 
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tition field strengths. A promising mechanism capable 
of driving magnetic helicity fluxes was introduced by 
Vishniac & Cho (|200ip , who found that in the presence 
of shear the flux follows isocontours of constant veloc- 
ity (Brandenburg & Subramanian 2005b; Subramanian & 
Brandenburg 2006 ) . Wc expect that this shear-driven mag- 
netic helicity flux is of crucial importance for large-scale 
dynamo action. 

The remainder of the paper is organised as follows: the 
numerical model is described in Sect. [2j and the results and 
conclusions are presented in Sects. [3] and respectively. 

2. The model 

Our model setup is similar to that used by Brandenburg 
et al. (|1996j) , Ossendrijver et al. (|2001[) and Kapyla et al. 
(|2004[ I2006|) . A rectangular portion of a star is modelled 
by a box situated at colatitude 6. The dimensions of the 
computational domain are (L x , L y , L z ) = (4,4, 2)d, where 
d is the depth of the convectively unstable layer, which is 
also used as the unit of length. The box is divided into three 
layers, an upper cooling layer, a convectively unstable layer, 
and a stable overshoot layer (see below). The following set 
of equations for compressible magnetohydrodynamics is be- 
ing solved: 
DA 
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where V/Vt = d/dt + (U + U ) ■ V, and U = (0, Sx, 0) 
is the imposed large-scale shear flow. A is the magnetic 
vector potential, B = V x A the magnetic field, and 
J = V x B / no is the current density, po is the magnetic 
permeability, r\ and v are the magnetic diffusivity and kine- 
matic viscosity, respectively, K is the heat conductivity, p 
the density, U the velocity, g = —gz the gravitational ac- 
celeration, and ft = £Iq(— sin 6*, 0, cos 6) the rotation vec- 
tor. The fluid obeys an ideal gas law p = pe{pf — 1), where 
p and e are the pressure and internal energy, respectively, 
and 7 = cp/cy = 5/3 is the ratio of specific heats at con- 
stant pressure and volume, respectively. The specific inter- 
nal energy per unit mass is related to the temperature via 
e = cyT. The rate of strain tensor S is given by 

Sy = i(C/ij + U j4 ) - |%V • U. (5) 

The last term of Eq. (j4]) describes cooling at the top of 
the domain. Here, t(z) is a cooling time which has a pro- 
file smoothly connecting the upper cooling layer and the 
convectively unstable layer below, where r — > oo. 

The positions of the bottom of the box, bottom 
and top of the convectively unstable layer, and the top 
of the box, respectively, are given by (z 1; z 2 , z 3 , 24) = 
(— 0.85, 0, 1, 1.15)rf. Initially the stratification is piecewise 
polytropic with polytropic indices (7711,7712,7713) = (3, 1, 1), 
which leads to a convectively unstable layer above a sta- 
ble layer at the bottom of the domain and an isothermal 
cooling layer at the top. All simulations with rotation use 
9 = 0° corresponding to the north pole. 



2.1. Nondimensional units and parameters 

Non-dimensional quantities are obtained by setting 

d = 9 = po = c P = po = 1 , (6) 

where po is the initial density at Zi. The units of length, 
time, velocity, density, entropy, and magnetic field are 



[x]=d, [t] = Vd/g, [U] 

[s] = c P , [73] = V dgpopo . 



% , [p]= po, 



(7) 



We define the fluid and magnetic Prandtl numbers and the 
Rayleigh number as 



„ v „ v „ qd 4 
Pr = — , Pm = - , Ra = - — 

Xo V vxq 



1 As 
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where xo = K/(p m cp) is the thermal diffusivity, and p m is 
the density in the middle of the unstable layer. The entropy 
gradient, measured in the middle of the convectively unsta- 
ble layer in the non-convecting hydrostatic state, is given 
by 



1 ds 

cp dz 



V- V 



ad 
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where V — V a d is the superadiabatic temperature gradient 
with V a d = 1 — 1/7, V = (cUnT/<91np) Zm , where z m = z 3 — 
Z2, and Hp being the pressure scale height (Brandenburg 
et al. I2005P . The amount of stratification is determined by 
the parameter £0 = (7 — l)eo /(gd), which is the pressure 
scale height at the top of the domain normalized by the 
depth of the unstable layer. We use in all cases £0 = 0.3, 
which results in a density contrast of about 23. We define 
the magnetic Reynolds number via 



Rm 



r/kf 



(10) 



where k{ = 2n/d is assumed as a reasonable estimate for 
the wavenumber of the energy-carrying eddies. Note that 
our definition of Rm is smaller than the usually adopted 
one by a factor 2tt. The amount of shear and rotation is 
quantified by 
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The denominators in Eq. (jlip give an estimate of the con- 
vective turnover time. The equipartition magnetic field is 
defined by 



B, 



eq 



(popuy/ 2 , 



(12) 



where the angular brackets denote volume averaging. 



2.2. Boundary conditions 

Stress-free boundary conditions are used for the velocity, 



u,. 



u z = o, 



(13) 



and either vertical field or perfect conductor conditions for 
the magnetic field, i.e. 

B x = B y = 0, (Vertical Field) (14) 
B x , z = B y . z = B z = 0, (Perfect Conductor) (15) 
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Table 1. Summary of the runs. The numbers are given for the saturated state of the dynamo (except for Run El, 
which was not run to saturation). Runs D and D3 are the same model. Here, Ma = u rms / '(gel) 1 / 2 , B lrns = B rms /B eq , 

(By) 1 / 2 = {W y y/ 2 /B Iins , and Q- 1 = (B^) 1 / 2 / (B^) 1 / 2 . The last column gives the magnetic field boundary condition at 
the z-boundaries, except for the last two lines where the x-boundary condition is given. Pm = 5 in all models except D6 
where Pm = 10. 
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respectively. We may think of them as open and closed 
boundaries, respectively, because they either do or they do 
not permit a magnetic helicity flux. In the y and x direc- 
tions we use periodic and shearing-periodic boundary con- 
ditions, respectively. The simulations were made with the 
Pencil CodeQ, which uses sixth-order explicit finite dif- 
ferences in space and third order accurate time stepping 
method. Resolutions of up to 512 3 mesh points were used. 

3. Results 

3.1. Description of the hydrodynamic state 

The hydrodynamics of local rotating convection simulations 
similar to ours has been discussed previously in the liter- 
ature (e.g. Brandenburg et al. 1996; Kapyla et al. 2004). 
However, the imposed shear flow used in the present runs 
is a new element. Thus it is interesting to consider the dif- 
ferences to the cases where it is absent. We consider four 
representative runs where physical ingredients are added in- 
dividually in order to see their respective effects. In Run A 
neither rotation nor shear are present, whereas in Runs B 
and C only rotation and only shear are applied, respectively. 
In Run D both shear and rotation are used (see Table Q] for 
details). 

All our runs start with a weak magnetic field (10 -5 
times the equipartition value). For Runs A to D it takes 
about 200 turnover times for this field to grow to values 
such that it affects the flow, so we consider this as the kine- 
matic stage of the simulation. We use this time span to 
characterize the basic hydrodynamic state. 

1 http : //www.nordita. org/sof tware/pencil-code/ 



The upper row of panels in Fig. [T] shows the vertical 
velocity in the middle of the convectively unstable layer for 
Runs A to D. In the cases without shear (Runs A and B), 
the convective pattern is dominated by large cells which are 
isotropic in the horizontal plane. In the case with rotation 
(Run B), the downflows at the vertices of convective cells 
tend to exhibit vortical structures, contributing to net he- 
licity. When shear is added the convective cells are slightly 
elongated along the y direction, which is the direction of 
the imposed large-scale flow. 

The horizontally averaged root mean square velocity 
shows a very similar profile in all models, i.e. peaking near 
the top of the convectively unstable layer and decreasing 
towards the base of the layer, see Fig. [5J The runs without 
shear are close to each other, with somewhat smaller veloc- 
ities in the case with rotation (Run B). When a horizontal 
shear flow is added, the rms velocity is somewhat larger, 
but in Run D u rms decreases again. The shear flow itself 
does not contribute to the definition of the rms velocity. 

The left panel of Fig. [3] shows the horizontal averages 
of kinetic helicity density for Runs A to D averaged over 
the initial kinematic stage of the simulations. In the case 
without rotation or shear (Run A) very little net helicity is 
produced. Adding rotation corresponding to the north pole 
(g ■ < 0; Run B) produces significant negative helicity 
within the convection zone, suggestive of a positive a-effect. 
In the case with only shear (S < 0, so that g ■ V x U > 0; 
Run C) the helicity profile is similar to Run B, but of oppo- 
site sign and roughly half the magnitude. If both rotation 
and shear are present, as in Run D, the sign and profile 
are very similar to Run B, but the magnitude is smaller by 
approximately a factor of two. This suggests that there is 
partial cancellation of the two effects. 
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Fig. 1. Vertical velocity in the middle of the convectively unstable layer of Runs A-D (from left to right) during the 
kinematic stage of the simulations [upper row of panels, t — 200 {d/g) 1 ^ 2 corresponding to tu Tms k{ « 50] and the saturated 
state [lower row, t = 2000 (d/g) 1 / 2 or tu rms k { w 500]. 




Fig. 2. Horizontal averages of the rms velocity in the kine- 
matic (left panel) and saturated states (right panel) for 
Runs A to D. Linestyles are indicated in the legend in the 
left panel. The vertical lines at z = (0, l)d denote the base 
and top of the convectively unstable layer, respectively. 

In the late stages of the simulations when the magnetic 
field has saturated the convective pattern does not change 
significantly in the cases without shear (see the lower row of 
Fig.[T]). In Runs C and D, there seems to be a larger change, 
i.e. clearer elongation of cells along the y direction, possibly 
due to the strong large-scale magnetic field that develops in 
those cases. Figures [5] and [3] show that velocity and helic- 
ity decrease somewhat in comparison to the hydrodynamic 
state in all cases except Run D where the helicity seems 
to increase. Part of this change can be attributed to the 
fact that some of the runs were not yet fully relaxed dur- 
ing the kinematic phase and that they are still undergoing 
a slow thermal adjustment, especially in the cases without 
shear. However, in Runs C and D where the dynamo grows 
to super-equipartition strengths, the magnetic field is likely 
to be responsible for much of the change. 

In the absence of shear the spectra are similar to those 
published previously (e.g. Brandenburg et al. |T996), but 



Fig. 3. Horizontal averages of kinetic helicity in the kine- 
matic (left panel) and saturated states (right panel) for 
Runs A to D. Linestyles are indicated in the legend in the 
left panel. 

with shear we should really only look at one-dimensional 
spectra in the streamwise direction, because this is the only 
periodic direction. In the kinematic stage little difference is 
seen between the four runs (see the left panel of Fig. 
However, in the saturated state the spectra for the shear- 
ing runs are steeper than for the non-shearing runs. This 
may indicate the development of large-scale ordered veloc- 
ity structures. 

3.2. Excitation of dynamo action 

We consider first the effects of shear and rotation on the 
dynamo. In Fig. [5] we show the resulting growth of the 
rms magnetic field for Rm = 56 ... 74 and Pm = 5 for 
Runs A to D introduced in the previous section. The slow- 
est growth occurs in Run A with no rotation and no shear 
(A = 0.016 u rms k{). Adding rotation (Run B; Co = 0.42) 
almost doubles the growth rate (A = 0.03 u rms fcf ), and the 
saturation field strength of the horizontally averaged field is 
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Fig. 4. Power spectra of velocity from the middle of the 
convectively unstable layer from the kinematic (left panel) 
and saturated states (right panel). Linestyles as in Fig. [5] 
The straight lines with a slope of —5/3 are given for com- 
parison. 
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Fig. 5. Growth of the total rms magnetic field from four 
runs, which were performed without shear or rotation (solid 
line, Run A), with only rotation (dashed, Run B), only 
shear (dash-dotted, Run C), and with shear and rotation 
(triple-dot-dashed, Run D), respectively. Pm = 5, Rm w 
60 . . . 90, and grid resolution 256 3 in all runs, see also Table 
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Fig. 6. Horizontally averaged magnetic fields B x (upper 
panel) and B y (lower panel), normalized by the volume 
average of B cq , as functions of time and z for Run D5 with 
Rm = 121, Sh = -0.21, and Co = 0.42. The dotted white 
lines show top (z = d) and bottom [z — 0) of the convection 
zone. 
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Fig. 7. Same as Fig. |6l but for Run C with Rm = 74, 
Sh = -0.17, and Co = 0. 



somewhat higher: B w 0.4 B cq . Next, adding shear (Run C) 
raises the growth rate further and, more importantly, it 
raises the saturation field strength by a factor of about 3 
(B « l.l_B cq ). Again, adding rotation (Run D) increases 
the growth rate further (A = 0.04 u lms kf ), and increases 
the saturation level by almost fifty per cent. In the absence 
of shear and rotation the critical Reynolds number is ~ 30 
for open (vertical field) boundary conditions (see Sect. 13.61 
for more details on the effects of boundary conditions). If 
shear is included, the dynamo is excited already at Rm « 5. 

We note that in the kinematic phase the mean and to- 
tal fields have the same growth rate. The two begin to de- 
part from each other only in the saturation phase which 
is achieved more quickly for the small-scale field. This be- 
haviour is reminiscent of dynamos where the saturation of 
the large-scale field is controlled by magnetic helicity evo- 
lution (Blackman & Brandenburg I200"2")) . 

3.3. Large-scale magnetic fields 

Adding a large-scale shear flow not only makes the dynamo 
easier to excite, it also helps the generation of large-scale 
magnetic fields and, most importantly, it is expected to 



drive magnetic helicity fluxes along lines of constant shear 
(Vishniac & Cho l200fl Subramanian & Brandenburg 2006). 
A representative solution from Run D5 with Rm = 121, 
Sh = -0.21, and Co = 0.42 is shown in Fig. El Note that 
in runs with rotation the large-scale field has opposite sign 
in the overshoot layer and the convection zone, whereas in 
models with only shear the field has the same sign every- 
where; see Fig. [7] 

In all cases the magnetic energy grows exponentially 
from a weak small-scale seed magnetic field. Generally the 
small-scale dynamo saturates first and a large-scale field 
develops later; see Fig. [8] for representative results from 
Run D5. Figure [9] shows that during this process the mag- 
netic field changes from a filamentary field to a more dif- 
fuse one. This behaviour has also been seen in simula- 
tions of forced turbulence with shear and open boundaries 
(Brandenburg 2005). The signs of B x and B y are always 
opposite to each other; see Fig. [TUJ This is simply because 
of negative shear, S < 0, so, for example, a positive B x 
results in a negative B y . In some simulations we have ob- 
served sign changes of the large-scale field in the saturated 
stage but these do not seem to be periodic but rather one- 
off events. 
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Fig. 9. Snapshots of B y in the early phase (left: t = 200 (d/g) 1 / 2 , corresponding to iit rms fcf « 50) and saturated 
phase (right: t — 1500 (d/g) 1 / 2 , or tu lms kf ss 360) for Run D5. The sides of the box show the periphery of the 
domain whereas the top and bottom slices show B y at vertical heights z = d and z = 0, respectively. See also 
http : //www.helsinki . f i/~kapyla/movies .html 
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Fig. 8. Rms- values of the total magnetic field (solid line), 
the mean field B lms = (B x + By) 1 ' 2 (dashed line) and 

fluctuating field with b 2 ms = B 2 ms — B lms (dot-dashed line) 
for Run D5. 

The rms value of the mean field contains a major 

fraction of the total field, i.e. (bI) 1 / 2 / B lms w 0.5... 0.8. 
Without shear the large-scale field is weaker by an order 
of magnitude (Table [T]). Note that at comparable Reynolds 
numbers the small-scale dynamo is excited in the absence of 
rotation and shear (Run A). A very similar large-scale field 
pattern is also obtained for Rm ss 10 and « 30 where the 
small-scale dynamo is absent or marginally excited. Further 
below we will show that the final saturation levels are only 
weakly dependent on Rm. 

It turns out that in the range 0.05 < |Sh| < 0.25 the 
growth rate is proportional to the shear rate, so the ratio 
A/ 1 5 1 w 0.1 is approximately independent of the value of 

\S\; see Fig.QTJ where A is the growth rate of (By) 1 / 2 . The 
same scaling was found in the simulations of non-helical 
turbulence with shear (Yousef et al. l2008aj) . It becomes in- 
creasingly difficult to perform runs with — Sh > 0.2 without 
rotation due to the simultaneous generation of large-scale 
vorticity. The theory for such a 'vorticity dynamo' has been 
developed by Elperin et al. ([20031) and has also been found 
numerically from forced turbulence simulations by Yousef 
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Fig. 10. Horizontal averages of the total rms magnetic 
field (top panel), B x (z) (middle), and B y (z) (bottom) for 
Runs A to D. The data are averaged in time over the satu- 
rated state of the dynamo. Linestyles as indicated in Fig.O 
The dotted vertical lines denote the convectively unstable 
layer between < z < d. 



et al. (|2008a[ I2008b|) and Brandenburg et al. ([2008]) . The 
vorticity generation is an interesting subject in itself and 
has been studied in a separate paper (Kapyla et al. [5008) . 
To avoid complications due to the vorticity dynamo we add 
rotation into the system which stabilizes the shear terms in 
the Navier-Stokes equations. 

When rotation is added, the growth rate divided by the 
shear rate seem to stay approximately constant for Sh > 
—0.25. When shear is increased further, — X/S decreases 
rapidly, see the left panel of Fig. [TTJ On the other hand, if 
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Fig. 11. Left panel: growth rate of (By) 1 / 2 divided by the 
shear rate S as a function of Sh for two set of runs with 
Co = (solid line) and Co « 0.5 (dashed line). Right panel: 

growth rate of (By) 1 / 2 normalized by the rotation rate f2 
for two set of runs with Sh « —0.25 (solid line) and Sh = 
— | Co (dashed line). Rm « 30 in all runs. The horizontal 
dotted lines show curves proportional to the shear rate S 
(left panel) and rotation rate fl (right panel) for reference. 

the shear is kept constant the quantity A/f2 stays constant 
up to Co w 0.5 after which it decreases rapidly, see the 
right panel of Fig. [TTJ The last case considered is to fix the 
ratio q = —S/n = 1.5. In that case the behaviour of X/il 
is almost identical to the case where Sh was kept fixed. A 
very similar trend was found by Yousef et al. ( 2008b)). 

The left hand panel of Fig. [T2] shows one dimensional 
power spectra of the magnetic field in the saturated state 
for Runs A to D. The runs with shear exhibit a steeper spec- 
trum in comparison to the cases without shear. However, 
it is clear, that the spectrum of the streamwise component 
of the magnetic field is sharply peaked at k y = and drops 
almost by an order of magnitude at the next wavenumber 
1 (in units of fci), because of the large scale field. This is 
the reason why we have also plotted spectra in linear scale 
in the right hand panel so as to accommodate the k y = 
wavenumber. 

3.4. Dependence on box size 

Convective overshooting in the present simulations is com- 
parable to that presented in Kapyla et al. Ij2004p and cer- 
tainly does not reach all the way down to the lower bound- 
ary. However, in order to determine whether the lower 
boundary still affects the results appreciably, we modify 
Run D2 such that the lower stable layer reaches down to 
z = — 1.85rf in comparison to z = — 0.85d in the other runs. 
We denote this run as D2-deep. We find that the large-scale 
magnetic field develops much the same as in the standard 
case, with the exception that the saturation of the field in 
the overshoot layer is slower, see Fig. [T3J In the saturated 
state the magnetic field occupies the whole overshoot layer. 
This is probably due to the large diffusion and the shear 
flow in the overshoot layer which are absent in the Sun. 

Figure [T] shows that relatively few convective cells are 
present in the computational domain. In order to test the 
effects of increasing the horizontal size of the domain we 
again take Run D2 as the basis and double the horizontal 
extent. This run is referred to as D2-wide. We find that the 
average magnetic field shows less fluctuations. This is due 
to the averaging over a larger number of convective cells. 
Growth rate and saturation level of the magnetic field are 



Fig. 12. One-dimensional power spectra of the magnetic 
field -Em(^) (left panel) and the y-component of the mag- 
netic field IwJ (k) (right panel) in the saturated state of 
the dynamo. Linestyles as indicated in Fig. [5] The straight 
lines with a slope of —5/3 are given for comparison. 
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Fig. 13. The rms value of B y for the runs indicated in the 
legend. The inset shows the same in linear scale. 

slightly larger than those of Run D2, see Fig. [131 This is 
because a bigger box results in more scale separation, which 
in turn increases the mean-field dynamo effect; see Eq. (80) 
of Brandenburg et al. (2002). 

3.5. Dependence on Rm 

Originally catastrophic quenching of the mean-field dy- 
namo effect was conjectured to derive from the action of the 
small-scale magnetic field whose magnitude was expected 
to increase as the magnetic Reynolds number is increased 
(Vainshtein & Cattaneo I1992[) . This should lead to Rm- 
dependent saturation levels of the magnetic field. In our 
Runs Dl to D6 listed in Table [IJ Rm is varied from 11 
(roughly a third of the critical value for exciting the small- 
scale dynamo) to 250 (roughly eight times supercritical). In 
all of these runs the saturation level of the magnetic field is 
essentially the same, bearing in mind the large fluctuations 
in the relatively short time series of the higher resolution 
runs; see Fig. [T3J Below, we demonstrate that the bound- 
ary conditions play a much more crucial role in allowing 
large-scale dynamos to operate efficiently. 

We recall that our definition of fcf in terms of the 
depth of the convectively unstable layer (2ir/d) is some- 
what arbitrary. However, it turns out that our values of 
A/(u rms fcf) are in agreement with those of forced turbu- 
lence (Haugen et al. 2004), where fcf is well defined. Over 
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Rm~ 11 (Run Dl) 
Rm~ 25 (Run D2) ' 
Rm~ 58 (Run D3) 
Rm~ 83 (Run D4) 
Rm~121 (Run D5) 
Rm~250 (Run D6) 
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Fig. 14. Upper panel: root mean square values of B y for the 
Runs Dl to D6. Linestyles as indicated in the legend. Lower 
panel: close-up of the interval where the dynamo saturates 
for the five largest Reynolds numbers. Data for the lower 
Reynolds numbers are shifted so that the saturation occurs 
approximately at the same point in the figure. The inset 
shows the saturated state in linear scale. 
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Fig. 15. Growth rate of the total magnetic field for runs 
with only shear (solid line) and Runs Dl to D6 with both 
shear and rotation (dashed line). 
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Fig. 16. Root mean square values of the total magnetic field 
and the horizontally averaged component B y for Runs C 
and C. VF denotes vertical field and PC perfect conductor 
boundary conditions. 



3.6. Effects of boundary conditions 

Here we compare Runs C and C made with vertical field 
and perfect conductor boundary conditions, respectively. 
The former allows flux of magnetic helicity across the 
boundary whereas the latter does not. Figure [TBI shows the 
time evolution of B rms and (By) 1 / 2 for these runs. The run 
with vertical field boundary conditions clearly exhibits dy- 
namically important large-scale magnetic fields. The rms 
value of the mean magnetic field is approximately 60 per 
cent of the equipartition value of the turbulent velocity field 
(see also Table UJ). Contrasting this with the same run per- 
formed with perfect conductor boundaries, the difference 
is striking: although the total magnetic field is roughly 40 
per cent of the equipartition value by virtue of small-scale 
dynamo action, the mean field is weak, contributing now 
only roughly 20 per cent of the total. 

The difference between the runs can be understood in 
terms of magnetic helicity evolution: allowing a non-zero 
vertical field at the boundary enables small-scale magnetic 
helicity to escape from the computational domain. With 
perfect conductor boundary conditions the magnetic helic- 
ity flux vanishes and thus the total magnetic helicity can 
change only resistively, which leads to slow saturation of the 
large-scale field. The fact that large-scale magnetic fields 
have not yet been seen in convection simulations with per- 
fectly conducting boundaries does therefore not mean that 
this is not possible, but rather that it would presumably 
take a very long time (see the slow but persistent rise for 
model C in Fig.EBJ). 



3. 7. Horizontal versus vertical shear 



the range 35 < Rm < 250, our values of A scale still almost 
linearly with Rm (see Fig. [TS"]) ; it is known from forced tur- 

bulence simulations that the Rm ' scaling occurs only for 
Rm > 200. 

The growth rate of the mean field in the late stages of 
growth does not seem to be strongly dependent on Rm, see 
the lower panel of Fig. [Mj Note also that the saturation 
field strengths are similar for different values of Rm. 



In an earlier study, Tobias et al. (2008)) used a vertical 

shear profile, i.e. E7 ° (z), and found that no appreciable 
large-scale magnetic fields were generated even though open 
boundary conditions were used in the z-dircction. In order 
to compare with their results we have made runs with a 
shear profile 



£7 (0) = $U 



1 + tanh f 



z - z 2 
di 



-2/ ' 



(16) 
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Fig. 17. Root mean square values of the horizontally av- 
eraged component B y for Runs VSA' (solid line) and VSA 
(dotted line). The inset shows the rms values of the to- 
tal magnetic field from the same runs. NF denotes normal 
field condition and P periodic boundary condition in the 
^-direction. 
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Fig. 18. Horizontally averaged magnetic fields B x (upper 
panel) and B y (lower panel), normalized by the volume 
average of B cq , as functions of time and z for Run VSA'. 

where U = -O.l(gd) 1 / 2 and d\ = QAd. This flow is im- 
posed through an additional relaxation term in the Navier- 
Stokes equation of the form 

d t u = ...-T- 1 {y-u m ), (17) 

where r = 0.5 (d/ g) 1 ^ 2 . In these runs with vertical shear we 
put S = 0. 

We confirm the results of Tobias et al. (j2008|) in the 
case where the horizontal boundaries are periodic and find 
that the rms value of the large-scale (horizontally averaged) 
field amounts to only a few per cent of the total field. We 
believe that this is due to the inefficiency of the magnetic 
hclicity flux through the vertical boundaries. In this case 
the isocontours of shear are horizontal and thus the shear- 
mediated flux, which we conjecture to be efficient, cannot 
leave the system. To test this hypothesis, we have made ad- 
ditional runs where the x-direction is no longer periodic but 
instead a stress-free boundary condition for velocity and a 
normal-field condition for the magnetic field according to 

U x = U y , x = U Z}X = B y = B z =0 (18) 
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Fig. 19. Vertical magnetic field component at the top of 
the convectively unstable layer at z — z% for Run D2- 
wide toward the end of the run. The plot is oriented 
such that the mean shear flow points to the left in the 
lower part of the plot (x/d = 4) and to the right 
in the upper part of the plot (x/d = —4). See also 
http : //www. helsinki . f i/~kapyla/movies .html 



is used. In this case the shear-mediated magnetic helicity 
flux can escape through the ^-boundaries and significant 
(up to 50 per cent of the total) large-scale magnetic fields 
are indeed produced, see Fig. [TTJ This result lends credence 
to the hypothesis that the shear-induced flux is an impor- 
tant ingredient in allowing large-scale dynamos to saturate 
to full strength on a dynamical timescale. 

In order to illuminate the circumstances around the 
time when the mean toroidal field becomes suddenly quite 
weak (tu rms k{ w 1600), we show in Fig. [T5] the mean fields 
in the cross-stream and streamwise directions. Note that 
there is no field reversal associated with the sudden drop 
of toroidal magnetic field, as one might have expected from 
simulations of dynamo action that is driven purely by mag- 
netic buoyancy, as in the work of Cline et al. (|2003j) . 

3.8. Surface appearance of large scale field 

In order to make contact with observations, one must even- 
tually compute the field as it would be observable at the 
surface of the domain. This process can be rather com- 
plicated and would involve radiation transfer of polarized 
light. Here we consider instead just the vertical components 
of the magnetic field at the top of the convectively unstable 
layer at z = z^; see Fig. [TH] 

Two things are immediately evident: firstly, the mag- 
netic structures tend to be inclined by up to 30 degrees 
relative to the y-direction and, secondly, there are bipolar 
regions with a systematic magnetic field orientation such 
that a negative polarity follows a positive one, which is 
in agreement with an overall positive mean toroidal field. 
Comparing with a similar plot of Brandenburg ( 2005), the 
bipolar regions are here much more compact. Also, the incli- 
nation is here such that it would correspond to the southern 
hemisphere, because dU v /dx is here negative. 
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4. Conclusions 

The present simulations provide a clear demonstration 
that convection can produce large-scale magnetic fields of 
equipartition strengths in an open Cartesian domain with 
shear crossing the surface. Our results also demonstrate 
that this is not possible within a comparable time span 
using a closed domain with perfectly conducting bound- 
aries. The sensitivity on boundary conditions has also been 
noted in earlier work with forced turbulence and shear 
(Brandenburg 12005]) . This can be interpreted as being due 
to magnetic helicity fluxes allowing the system to dis- 
pose of excess small-scale magnetic helicity along contours 
of constant shear (Vishniac & Cho 2001; Brandenburg & 
Subramanian 2005b; Subramanian & Brandenburg 2006). 
In our case this is the vertical direction, but in Tobias et 
al. ([5008), it is the horizontal direction. Thus, in our case 
shear-driven magnetic helicity flux can escape the domain 
whereas in theirs it cannot due to the periodicity in the 
horizontal directions. This could explain why no large-scale 
dynamo develops in their case. Their model is otherwise 
similar to ours (same values of Co and Pm, whilst Rm is 
about two times smaller, and Sh is twice as large.) Control 
simulations show that replacing their periodic boundary 
condition in the x direction by a normal-field boundary 
condition allows strong large-scale fields also in their setup 

(Fig. En). 

The precise nature of the dynamo found here cannot 
be pinned down rigorously unless one is able to identify 
unambiguously the mechanisms that are actively at work. 
In addition to the er-cffect, there are other possible alter- 
natives, in particular the incoherent alpha-shear dynamo 
(Vishniac & Brandenburg |1997l Proctor r2007p . which could 
produce a large-scale magnetic field even without stratifi- 
cation and hence no regular a-effect. This interpretation 
has been favoured in some recent work where the shear- 
current effect was found not to be excited (Brandenburg 
et al. l2008j) . However, this matter is still under debate and 
there are arguments in favour of a shear-current effect (e.g. 
Kleeorin & Rogachevskii l2008p . 

In order to disentangle the relative roles of different 
mechanisms in the present work it will be important to 
establish a closer connection with mean-field theory by de- 
termining the components of the a and turbulent magnetic 
diffusivity tensors for different boundary conditions and as 
a function of Rm. Only when the results are shown to be 
converged with Rm may we expect them to be relevant for 
the Sun and perhaps other astrophysical bodies. 

We emphasize that large-scale dynamo action is nor- 
mally only seen in the nonlinear stage, whilst the linear 
stage is dominated by small-scale magnetic fields. It is 
therefore important to consider mean-field transport co- 
efficients that are affected by the magnetic field. It is in 
principle even possible that the relevant dynamo mecha- 
nism is an intrinsically nonlinear one, as it is in the case 
of accretion disc turbulence where the turbulence is the re- 
sult of the dynamo itself (Brandenburg et al. I1995[) . The 
present results do not give any indications in this direc- 
tion, because there seems to be a continuous transition from 
low-Rm to higher-Rm large scale dynamo action. For ex- 
ample at Rm =11, small-scale dynamo action in the usual 
sense is not expected, and nevertheless, the final satura- 
tion strength of the large-scale field is similar to the case of 
highcr-Rm dynamos. 



Another important aspect to consider in future work 
is the degree of scale separation. As the Reynolds number 
increases, the turbulence becomes more vigorous and the 
convection plumes will no longer penetrate the entire depth 
of the unstable layer. Therefore, we expect the wavenumber 
of the energy-carrying eddies to increase, i.e. kf > 2ir/d. 
This means that the degree of scale separation, i.e. the ratio 
fcf/fci, will also increase. However, as discussed at the end of 
Sect. 13.41 more scale separation is only expected to enhance 
large-scale dynamo action. 
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